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ABSTRACT 


An examination and a comparison of the relative merits of 
the finite element and Fourier series methods of solving radially 
loaded circular ring problems are made. The procedure employed 
to evaluate the two methods is to use each method to solve for three 
different load conditions and to compare the performance of the two 
methods on the basis of accuracy, ease of usage, and equipment 
required. The results indicate a айй солу accuracy for Both 
methods under most conditions. The Fourier series method is 


superior for solving problems with a distributed load condition. 





The finite element method is superior for solving problems with 


concentrated loads. 
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США TER I 


INTRODUCTION 


The structural analysis of many aerospace vehicles, suchas 
rockets, re-entry bodies, aircraft, space capsules, etc., is 
frequently accomplished using either the Fourier series method 
or the finite element method in conjunction witha digital com- 
puter. A considerable amount of effort has been devoted to the 
development and refinement of these two methods, but little 
attention has been given to a direct comparison of the methods. 
Fundamental questions regarding the similarities, differences, 
relative accuracy, ease of usage, and advantages or disadvantages 
of each of the methods have not been answered as yet. The ob- 
jective of this study was to seek answers to some of these funda- 
mental questions. 

The procedure employed to accomplish the objective was to 
apply the two methods of analysis to a portion of a typical aero- 
Space structure and to investigate the relative merits of the two 
methods. The portion selected for analysis was a thin circular 
rie with a rectangular cross section. The ring was selected 
because it represents, ina simplified manner, the circular 
cylinders that are frequently found in aerospace vehicles. Three 
different load conditions were applied to the ring for each of the 
two methods of analysis so that the effects of the load distri- 
butions on the performance of the two methods could be investi- 


gated. The extent to which the analyses were carried out in 


15 


terms of numbers of significant figures carried, nurnbers of 
terms evaluated in the Fourier series method, number of 
elements used, and use of double precision techniques in the 
FORTRAN computer programs, was greater than that normally 
used in common engineering practice. The reason for such 
depth was to bring out any subtleties peculiar to either of the 
two methods that might not appear if normal engineering accu- 


racy Wem@e used. 


The author gratefully acknowledges the guidance and 
assistance of Professor Robert E. Ball in the preparation 


of this thesis. 
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CHA PTER П 
DERIVATION OF THE GOVERNING EQUA TIONS 
FOR THE FOURIER SERIES ANALYSIS 

Consider the thin ring shown in Figure l where a is the radius 
of the ring middle surface; b, the ring width; and t, the thickness. 
The angle ( is the generator of the ring and is positive in a counter 
clockwise direction. All points in the ring are located by the 
coordinates Q and z, where z originates at the middle surface and 


is positive in the outward radial direction. The normal force Ng? 


S Nç 
Not 55 “ὁ ia 3% 


З. 
Ээ каа За as 
/0 A 


N 


Ф 


à. LZ 


FIGURE 1 


c / 


RING CONFIGURA TION 
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the shear force Qo: the bending moment My: and the applied 
radial loading P, are shown in the positive direction. The dis- 
placements, also shown in the positive direction in Figure 1, are 
w in the radial direction and v in the tangential direction. 
Assuming that the ring retains its circular shape after de - 


formation, the equations of equilibrium can be written as 


$ М 
а — штэ О (2-18) 
d © 


at, — № = 29% (2-18) | 


ў Ф | 


| 
O 


ата 


iN, 


-Q, = O (2-1c) 
40 


Using Equation (2-1a) to eliminate Qo from Equations (2-1b) and 


(2-1c) leads to the two equilibrium equations 
d Na d M, 
a— - — = 
i$ 4 (2-24) 





(νι 


The circumferential strain, & Ф” at any distance z from 


the middle surface, is given in Refemence l as 


dv 2" 1 496 wW 


-a мыш с — — — amn w som 


+ 
9 ado a (a+z) yo? a +Z 








(2-3) 


This relationship is based upon the assumption that points on a 


normal to the middle surface prior to deformation remain on the 
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normal after deformation and that the thickness of the ring 
remains constant. Substituting Equation (2-3) into the plane 


stress form of the elastic constitutive law leads to 














РҮ: iv 221 4ш UJ 1 
Цэн add a (a+z) ¿q? а + 2. шан 
where σφ is the circumferential stress at the location (0, 2.) 
and E is the modulus of elasticity. 
The normal force and bending moment are given by 
t 
3 
№ 3 To dz (2-5a) 
Ps | 
2. 
and 
t 
2 
Me а To zdz (2-5b) 
E 
2 
Substituting Equation (2-4) into Equations (2-5a) and (2-5b), 
integrating and applying the limits yield 
7. 
V wW 
bc ο ο ыа ц š + W (2-6a) 
ф 3 3 
« [ d$ al dö 
and 
42 
М, = w t (2-6b) 
Φ a" dp” : 


where D- Et and K-Et?/12, ж 


*In order to use this form of D and K the natural logarithmic 
terms that result from the integration of Equations (2-5a) and 
(2-5b) must be expanded in a series in terms of t/2a where terms 
to the fifth power and above are dropped. 
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The problem has now been reduced to one with four unknuwns, 
A М №, and Mp, and four equations; two of the equations are 


obtained from the conditions of equilibrium 


le cow 
“аф T 


-0 (2-2a) 


4^M 
2 ф _ й 
аМь-а γέ =»; (2-2Ь) 


and two are obtained from the application of the elastic consti- 





tutive law and the strain-displacement relationship 


-pU К [аш (2-68) 
No= Eu et "up r 
d'w 
Mo= "ar E dier 3 Ы ы (2-6b) 


Taking the appropriate derivatives with respect to 0 of Equations 
(2-6a) and (2-6b) and substituting the results into Equations (2-2а) 


and (2-2b) yield 





<1 dv + w = © (25725 
ddl do 
and 
29 ω | atf, = 
E + Ú) ελ ων 28 | sO (2-7b) 


where k — ё/12а7. Equations (2-7a) and (2-7b) are the governing 
differential equations for a thin ring under the assumptions that 


have been made. 
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CHAPTER III 
DEV ELOPMENT OF THE FOURIER SERIES SOLUTION 
TO THE GOVERNING EQUATIONS 
The starting point in the development of the Fourier series 
solution to Equations (2-7a) and (2-7b) is to assume that the 
applied load, P,, is a general function of 0, symmetrically dis- 
tributed about two perpendicular planes, one through ( — 0 and 
0- ff, and one through 0-117/2. Thus, the load can be repre- 


La 


sented by the Fourier cosine series * 


со 
P [^ COS VA Фф (3-1) 
m=0,2,4 


where the coefficients P, are given by 


fT 
= н] Р, 46 (3-2а) 
d 


an 


- | 
= 31 4 Cos Modo (3-2b) 
O 


The radial displacement, w, which is an even function of ў 
for the ring under the load distribution considered, is also written 


in a Fourier cosine series as 


CO 
B=) Wy, СОЗ Y (3-3) 
Mm = 0, 2,4 





* The symmetry conditions on P, lead to the elimination 


of all odd values of m, i.e., P,,= 0 for all odd values of m. 


РЭ 


The tangential displacement v is an odd function of and is 


given by 


ОО 
Varia Vara Saa үү\ ф (3-4) 
м = 2,5 


The governing differential equations are used to obtain WEE 
and v, as a function ої Р. For the axisymmetric mode, 
where m= 0, the substitution of Equations (3-1) and (3-3) into 
Equation (2-7b) yields 

лаг 


О (1+ А) 


6 = (3-5) 


For the general case of m —0, substitution of Equations (3-3) 
and (3-4) into Equation (2-7a) results in 
2. — - 
W; V. a m Lge = О (3-6) 


Substituting Equations (3-1), (3-3), and (3-6) into Equation (2-7b) 


gives 
2. 
Wm = ---ᾱ--- 
MD (2-1) 
Therefore, according to Equation (3-6) 
р А р 
а 
Vm = - 


Х D vv (wm? - 1)* 
Substituting these expressions for w,, and vg, into Equations (3-3) 


and (3-4) leads to. 


а? Р, à 9 Ри соз wd 


— A mE н (3-7a) 
ра b 9 (νη; - 1)” 


J 
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в“ -- = 


\ = — ə ə—O  .n м 
A^ v ^ S 2 тү” 


(3- 7b) 


All other quantities of interest are determined by using these 
expressions for w and v in the appropriate equations. For 


example, according to Equations (2-6a) and (2-6b), the normal 


force is given by 


22 P COS Ул Ф 
No = at, = а. У. — ye — (3-8) 
W = L, t 


And the bending moment is 


2. ў, Р co Р соз мф 
М,- 1+} a 2. mi Ri (3-9) 


M = 2,4 


To permit a direct comparison between the Fourier 
series and finite element analyses of the ring, w, №, апа My. 
are selected for evaluation by both methods at the points ф = 0, 
30, 60, and 90 degrees. In Chapter IV, Pa and the specific 
equations for w, No: and M, are developed and evaluated at 


these values of ў for three different load conditions. 
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С НАР ТЕНЕ ГҮ 


APPLICA TION AND RESULTS OF THE FOURIER SERIES SOLUTION 


The three loads selected for use in this study are shown in 
Figure 2. All three loads are radial loads, directed inward, and 
are symmetrically distributed about two planes, one through 0 — 0 
and ф = TT and one through 0 — 1 1/2. Caselis a "line" or "knife 
edge" load of P pounds distributed over the width of the ring at 0= 0 


and ( = т. 





P/2«abí(lbs/qm) 
P/2 a b (lbs/m*) 
Case 1 «= 0 degrees 
Case 2 e 2 30 degrees 
Case 3 “= 60 degrees 
FIGURE 2 


LOADS CONSIDERED 


^2 





Case 2 and Case 3 loads are uniform pressure loads, each with a 
total value of P pounds covering the areas shown in Figure 2. The 
uniform loads are initially represented as P pounds distributed 


over the angle 2™ so that the surface load P, is given by 


У ТАУ __ ай 
в 2с«0а0 дый, 


Thus, substituting Equation (4-1) into Equation (3-1) gives 


Р 


a ee "4 


W = ο, ἡ 


Accordingly Equation (3-2a) becomes 


P Ё Ч 
ᾱ-- abl | “1. Ἡ 
O T- 


- — (4-2) 
ο Tab 


Similiarly, substituting Eguation (4-1) into Eguation (3-2b), inte- 


grating and applying the limits yield 


P 


P =- ————_[Sin m& — sin m(T—&)| 
m ff alo mex 
Therefore 
= 2 P sin mec (4-3) 
m Tab me 


for m even. For odd values of m 


= 


m 


as originally noted in Equation (3-1). 


23 


Applying Equation (4-2) and (4-3) to Equations (3-7a), (3-8) 


and (3-9) leads to the Fourier series solution for w, N. and M 








9 
р OD 
T Wr 5 соз эф (4-4а) 
ff b D (i) MM (m?- 
P ре sinmax COS мф 
т (4-45) 
Ф Γον ть ec Y 1) 
m= 2,4 
and 
Mp- B а АР + 2aP у Si me cos mó (4-4c) 
ть4+4) Th e< yñ ( vn 2 - 1) 


m = 2,# 
The non-dimensional forms of these three equations, denoted by 
W, No and Me, and obtained by using 11 t^/12a*, D = Tš asa d 


Е = b are 











. αυ. l 22 sinmx cosmo 
112 š = — ия = T = теш (4-5a) 
aM fT a? A (14 4) 27 - 

- ЫН ji © бүүр үүлс COS ^ 

М = a = — <= E ў A (4-5b) 
Р mz 2,4 і 

апа 

ni eMe "n 2 €, sinme cos m¢ 

MEC aea a чк = A ч ае; 
a P Та? А(4+ А) ft e mM (v? - 1) 


- 6 4 


24 


These expressions are applicable only to Case 2 and 3 loads. 
For Case 1 load, cx — 0. Thus, taking the limit of Equations 


(4-5a), (4-5b) and (4-5c)as cx —=0 leads to the solution for 





Case 1 load 
- di 2 ® cos MO 
w= = e İy 21  (m*-1)* (4-6a) 
Ta“ A(1+4) anes a 
— oe aue caa ш (4-6b) 
Ф T ТГ m? - 1 
and M= 2,4 
EC 2, 8 
Мь= - і - 5.4 cos mo — 
(a^ Alith) T m? - 1 


WA = Let 

The term - 1/1 a^ A (1+ k) that appears in w and Mç w111 now 
be examined for its contribution and significance. The expressions 

M- N*a (dO) ә ы U,- Malad) *a (do) 

| 2 АЕ гал” 

give the axial strain energy and the bending strain energy respec- 
tively of an element b ta(d$) of the ring. They show that the axial 
Strain energy is proportional to 1/A and the bending strain energy 


is proportional to 1/I. Using U 7 1/A and U. = 1/I leads to 


b 
U,/Uy z I/A. In Reference 2, the axial strain energy for such a 


curved element is shown to be very small compared to the 
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bending strain energy of the element provided that t/a is small. * 
Thus, the term - 1/T a^ A (1*k) in w and м, is related to the 
extension of the middle surface of the ring and should be small, 
If the middle surface is considered inextensional, the axial 
Strain energy of am element of the ring is zero»sand tbe term 

-1/ Т ад A(l*k) which is proportional to U, /U,, does not 
appear in the expressions for w and Ms. 


the magnitude of the contribution of the axial strain energy to the 


In order to determine 


solution, a t/a ratio must be selected. As an example, a ratio 
of t/a-z 1/10 is chosen and w and M, are calculated with and 
without the term containing I/A at ф= 0, 30, 60 апа 90 degrees 
for all three load conditions. The results of these calculations 
are found in Table I, where the values referred to as "EXT" 
include the -1/f a“ A (1*k) term and the "NON-EXT'' values do 
not. As is evident from Table I the differences between the "EXT'' 
and "NON-EXT'" values are very small, or nonexistent for the 
accuracy retained, when t/a = 10. 

An important consideration in the application of the Fourier 
series method is the number of terms of the series required 
to give a desired accuracy in the solution, i.e., how fast does 


ale 
” + 


This statement assumes there are no distributed loads on 
the Clement, 1.e., only g@oncenmtravccedeads ane consaaered,  Phis 
thesis considers distributed loads, and thus the statement is not 
completely applicable. 
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DABLE I 


FOURIER SERIES EXTENSIONAL AND NON-EXTENSIONAL RESULTS 


ἘΝ CASE 1 LOAD CASE 2 LOAD CASE 3 LOAD 


—- = —— ρα E Wo 











ç s Ext. Non-ext. Ext. Non-ext 
| -0.0744|-0. 0599 | -0.0596 | -0. 0289 -ο. 0281! 
| 0.0 -0,1289 | =0:128- | онә (10239 
| 0.318 | 0.190 0. 190 0. 079 0. 080 
0152 |-0. 0149 
30° | | 271 amd 
| 0. | ‚ 047 | 0. 048 
w | 0.0361 x 0. 0364 | 0.0295 | 0.0298 | 0.0141 | 0. 0143 
60” | М -0. 433 ‚-0. 433 |-0. 413 -0. 413 0.358 0. 358 
м, -0. 115 |-0. 115 aps Ies lowo ои 
w | 0.0680 | 0.0683| 0.0571 | 0.0574 | 0.0295 0. 0298 
90? No -0.500 |-0.500 |-0.477 |-0.477 1-0.413  .|-0.413 
! Mọ! -0.182  |-0.182 -0.159 | -0. 159 | -0. 095 -0. 095 
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the series converge? This $üformation for the inextensional 

form of the solution is shown in Tables II, Ill and IV for the first | 
4 terms of w, Мапа My at d= 0, 30, 60 and 90 degrees for all | 

three load conditions. The number of terms required for con- 

vergence of a series to four decimal place accuracy for the 

displacement and three decimal place accuracy for the forces 

and bending moments is shown in the last column with the con- 


verged value of the series. 
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TABLE TI 


FOURIER SERIES NON-EXTENSIONA L CONVERGENCE 
FOR CASE 2 LOAD 


4 Г Сопуегреа | 





Value.  , 
‚0744 ағ о = 14 | 
. 000 atm = 53% 
. 346 арт = 894 
. 0334 at m = 6 
. 230 etmez se 
.068 atm = 52 
0364 at“m = 14 
. 433 atm = 24 
. 115 atm = 38 
0083 ALT 8 
-0 500 at rm = 26 


-0.182 atm - 38 
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TABLE IH 


FOURIER SERIES NON-EXTENSIONA L CONVERGENCE 
FOR CASE 2 LOAD 


Converged 








Value 


.0596 at m 





125 atm 





.190 at m = 20 
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TABLE*IV 


FOURIER SERIES NON-EXT ENSIONA L CONVERGENCE 
FOR CASE 3 LOAD 





oo ши - e ——' ч а «νο. хс сше сш чч - as, 


Converged 


Value 





— . 


! 
e 
e 
N 
Qo 
-J 
b m 
Gh 4 
! 
нь 







-0.231 - 0. 239 


y 
| 
t 
i 
! 


|! 
pe 
e 


-0.239 atm = 





0. 088 0. 079 1 0. 079 : 0. 080 0.080 at m = 14) 
мете н 
HR ZUR 102786  1-0.278 Faen КЕТ Sai 33 

EE ums Ίσως [ams T o Oai 
0.0146 0. 0143 | 0. 0143 at m= 4 
-0.362 |-0.358 Ф 5 т 4 

: | 

-0. 044 |-0.039 |-0.039  |-0.040  |-0.040 atm= 8 
0.0292 | 0.0298 | ` 0.0298 atm= 4 
-0.406 |-0.415 |-0.415 |-0.414  '|-0.413 аб т 44 
Мор |-0. 088 |-0.097 (1-0.097 |-0.095 · |-0. — 8! 
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CHAPTER V 


DEVELOPMENT OF THE FINITE ELEMENT SOLUTION 


The finite element method is based upon the concept of 
replacing the actual continuous structure with a number of struc- 
tural elements of finite size, leading to an assembled structure. 
The oe сун of equilibrium and continuity at the element 
joints or nodes lead to a set of simultaneous algebraic equations. 
There are two general methods for formulating these equations, 
the displacement method and the force method. The displace- 
ment method, also known as the stiffness or direct stiffness 
method, is the one used in this thesis to solve the ring problem. 
This method is discussed in detail in Reference 3. A brief ex- 
planation of the method is given here for continuity. 

In the displacement method the displacements and rotations 
at the nodes of the elements comprising the structure are taken 
as the unknown quantities. The matrix equation which See 
the external forces to the tape elemi of an assembled structure 


is 
(Xi ік” lut в {хо} (5-1) 


Matrix ixi is the column matrix of nodal forces, i.e., forces 
that act upon the nodes of the element. Ina general three dimen- 
sional structure these may include bending and twisting moments 
as well as rectilinear forces. The matrix [2 is a square 
symmetric matrix for the entire structure. This matrix is 

known as the stiffness matrix and is assembled from the individual 


stiffness matrices of the elements forming the structure. The 
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entries in the element stiffness matrices are called stiffness 
influence coefficients. They are determined by displacing a 

node of the element by one unit in the direction of one of. the 
degrees of freedom permitted for that node and calculating 

the nodal forces required to maintain that unit displacement and 
to prevent movement in any other degree of freedom for the 
nodes on the element. The amount of force required at a node 

to maintain this displacement state is one stiffness influence 
coefficient. The complete set of stiffness influence coefficients 
is determined by repeating this procedure for each degree of 
freedom at every node. The column matrix (lu) is the matrix of 
displacements and rotations which may take place at a. node. The 
column matrix (x?) is a matrix of initial forces. They are deter- 
mined by clamping all of the nodes of the assembled structure 
and calculating the external force required at each node to 
maintain this zero displacement state under the applied load con- 
dition. Thus, loads distributed between nodes and concentrated 
loads not located at nodes can be treated. 

There are various techniques available for the manipulation 
and solution of Equation (5-1). The one used in this study is:to 
transpose (x?| to the left-hand side of Equation (5-1). This leads 
to a single column matrix of forces | {x} 5 Ра) ^ Meet, 
the matrix equation is reduced in order by a consideration of the 
boundary conditions on the structure. If the boundary conditions 
specify zero displacement at some nodes, then those nodal dis- 
placements are no longer unknowns and are removed from (а) 


along with the row of the equation and the column of [к] that 
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correspond to that zero displacement. The remaining matrix 
equation, which is referred to as the reduced matrix equation and 


is subscripted with (^, is 
ха - xo = [Ke] fur] (5-2) 


The unknowns in Equation (5-2) are contained in the matrix of 
displacements up) . Equation (5-2) is à set of simultaneous 
linear algebraic equations that is solved on a digital computer 
using the FORTRAN IV subroutine ''DSIMQ," a description of 
which appears in Appendix C. The outputs of 'DSIMQ'" are the 
nodal displacements {up} ' 

In order to determine the nodal forces [x| , the stiffness 
matrix equation is solved for each element. These equations, 


subscripted Ч e 


O ` 


ха - СО (5-3) 


` 


The unknowns in this equation are in EA since the displacements 
Í un) are obtained from the appropriate locations in m and (uj. 
The matrix (Хү) is obtained by adding the matrix (xil , Which 
is determined by selecting values from appropriate locations in 
О ў r P y ' . 

х } , to the matrix EN UT . With (un! and [X4 known, 
all of the nodal displacements and forces have been determined, 
and the problem is solved since the forces and displacements at 


any interior location of an element can be determined from the 


applied load and nodal forces and displacements. 
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Applying the linite element method to the ring problem under 
consideration requires the selection of the element geometry. 
Either straight or curved elements can be used. Curved elements 
are selected since they more naturally match the shape of the ring 
and result in displacements and forces that are readily interpreted 
without requiring a change in coordinates. The curved element 
used is drawn in Figure 3 with the nodal forces Му Qu N. Q. 
bending moments Μι, M5, displacements Жүр MET W У» апа 


rotations @,, 0, shown in the positive sense. The element covers 


à 
the angle 8 . A stiffness matrix for this element has been 


derived in Reference 3% based on the assumption that the thick- 


ness of the element, t, is small compared to the radius of 


Naa V 
di Rum 


Καν ΜΙ; 9 Θ; 


Qu^ 
яг 


м, Ө, 


а 
DN 


PIGURE 3 


FINITE ELEMENT CONFIGURA TION 





*А more detailed derivation and a comprehensive study of 
the curved element appear in Reference 4. 


5b 


curvature, a, so that only the strain energy of bending need be 
considered during the derivation. This assumption neglects the 
shear distortion energy, the energy resulting from any coupling 
between bending moment and normal force, the energy of axial 
deformation, and the displacements that correspond to these 
three forms of strain energy. The stiffness matrix equation from 
Reference 3 for the curved element can be given in a form corre- 


sponding to Equation (5-2)as 





/ О 
М, ЫГ М, 
ο 
Q1 “чі | ΩΙ 
M,/a a0, | Μι /8 
ED [y и (5448) 
N рац | N У Ne 
2 a 2 2 
ο 
92 = 22 
М,/а | 20, | M5* /a 


where [Ky Jis the six by six matrix of stiffness influence 
coefficients for the curved element. The matrix [ки |, each 
term of which is a function only of 3 ; willijnot be written out i 
functional form because it is quite lengthy. For this thesisa 
double precision digital computer subroutine, 'RELM,'" was 
written which computes | Ку | ох an input of 8 . A copy of 
"RELM" is included in Appendix C. The stiffness matrix EN 
con pured by" hPa" for f? — 7.5, 15, апа 30*"degrees-1is 
presented in Tables V, VI and VII respectively. 


With ES determined, ES can be assembled. Equation (5-1) 


is completed for the ring by using the loading conditions to obtain 


{κ} and the appropriate portions of {xl . Because of the 
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TABLE V 
к] for [3 — 1.5 degnees 
MIOS 12559962 26,675 -18, 65764; 1, 20860 
|, 223,063 85/518 2099 - 1522080002, 74, 855 
26,675 2099 69 2161, 2 | 1,400 


1816657520 -1, 223, 762 -26, 721 18, 665,685 -1, 225,063 





1, 225, 762 74,855 1,400 -1,223, 063 85,518 
| 
габ, 721 -1, 400 -23 267675 зе, 1099 69 , 
— — 
TABLE VI 
! 1 for ВЕ 15 deguees 
576,897 75, 862 3,906 -576,874 76,035 
MB , 662 10,658 524 -76, 035 95 339 
3, 306 524 34 -3, 329 340 
-576, 874 -T035  -3,829 516,099" -15,862 
76, 035 9, 399 350 -75,862 . 10,658 
-$, 52 -350 211! 3, 8/96 -524 
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TABLE VII 


EN а f? — 30 degrees 





17248 4,598 TRA 4,640 "121 
4, 598 ipie 7 130 -4,640 " 158 87 
399 130 17 -411 87 -6 
-1* 828 -4, 640 -411 17, 243 -4, 598 399 
4, 640 1, 158 87 -4, 598 1,887 EET 
| -an -87 -6 399 NB 0 17 


515) 





symmetrical nature of the applied loads only one quadrant of the 
ring needs to be considered. The quadrant from ў -"0 to 999 degrees 
is selected so that a direct comparison with the Fourier series 
solution for w, Ny and м, can be die. for ф= E 60 аша | 


90 degrees. 
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CHAPTER VI 


APPLICA TION AND RESULTS OF THE FINITE ELEMENT SOLUTION 


In order to apply the finite element solution to the ring 
quadrant Ó = 0 to 90 degrees, Equation (5-1) is assembled for the 
entire quadrant using Equation (5-4) for each of the elements com- 
prising the quadrant. [Ky] matrices for [3 = 7.5, 15 and 30 de- 
grees are used in this thesis so that the quadrant is divided into 
12, 6 and 3 elements respectively for the solution of each of the 
three load conditions. This assembled equation is reduced and 
solved in the manner described in Chapter V. The process of 
assembling, reducing, and solving Equation (5-4) is best described 
by uSing an example. Case 2 load with З — 30 degrees is selected 
as the example. 

Figure 4 shows the assembled quadrant and the node numbering 


sequence for the three elements in the quadrant. For three 


A 


FIGURE 4 


3 ELEMENT QUADRANT FOR CASE 2 LOAD 
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elements κ] is a 12 x 12 square symmetric matrix which is 
assembled using three of the 6 x 6 EN matrices from Table VII. 


When assembled for this example, Equation (5-1) appears as 








М) Ve 0 NO 
Q), = 0 | -У ] $i 
Mıla | LN o M, /a 
N,-0 | 2 x М)” 
| Bay Ὁ m pa | 
M,/a = 0 «ΕΙ ΚΙ 12x12 / a8, k. Qa llaq (6-1) 
N3 = 0 | V3 | ΤᾺ | 
IM FP | "3 1937 
M, /a ад, M,°/a 
Na Уд =0 Na” 
w | my 0.) Ол” 
Lal μα: ; αὶ αθ, Ξ 0 İM, /a 
| 4 ) \ 4 ) - 


where the zero entries in the (xi and fu] matrices are 
obtained from the symmetry characteristics of the load. The 
initial forces in ү are derived in Appendix B and presented in 
Table BI. They are transposed to the left-hand side of Equation 


(6-1) to give 
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N} - . 12709678 





0 <4 44822182125 1 
М,/а -  .01944670 
0 -  .12707618 
0 -  .22188951 
0 -  .019446710 
Л: = (6-2) 
0 
0 
Ял 
0 
| M,/a 


The reduced matrix equation is created from Equation (6-2) by 
striking out the rows of the equation and columns of [к] 
12x 12. The reduced equation with [Кр an 8 x 8 matrix 
appears in Table VIII. 

A computer program, "ВЕТА 30," a copy of which appears 
in Appendix C, solves Equation (6-3) for (up) by using the 
subroutine "DSIMOQO.'' For this example the output of !'DSIMO" is 


== ‚ № 
Wg == М е” 1 ac - 0. 0596 


DET T 0. 0255 
== -0. 0287 


ο Е мава? 
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TABLE IX 


FINITE ELEMENT RESULTS FOR CASE 1 LOAD 
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TABLE X 


FINITE ELEMENT RESULTS FOR CASE 2 LOAD 


D Lee Д аз | 
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& = 8 


-0. 0596 
-@ 126 
0.190 
“а 0287 
-0 2397 
0. 080 
0.2918 
-0. 413 
- 0. 095 
-0. 0573 
-0.477 


*0. 159 


TABLE XI 


FINITE ELEMENT RESULTS FOR CASE 3 LOAD 


ЭХ: 8 =з 





-0. 0287 -0. 0287 -0. 0287 
-0. 239 | анд | -0.239 
080 _ 0.080 0.080. 
0149 С -0. 0149 | -0. 0149 
i 30897] | 20. 2711 
.048 — ' 0.048 | 0.048 
“таз - ‘| noms | 0.0143 ' 
.358 858 -0. 358 
040 -0. 040 -0. 040 
к ‚0298 | 0. 0298 -0. 0298 
90°. «0.413 | -Q. 413 -0. 413 
м | -0095  '  ' -0.05 | = [ 
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CHAPTER VII 


COMPARISON AND DISCUSSION OF THE RESULTS 


In order to provide a standard to which the results of the 
Fourier series and finite element methods can be compared, the 


ring problem is solved analytically for N, and м, in Appendix A 


р 
for the three load conditions. The results of this analytical 
solution are non-dimensionalized and are entered in Tables XII, 
XIII and XIV with the results from the Fourier series and finite 
element methods. The Fourier series values in these three 


tables are the converged values of w, N, and Му, while the finite 


Ф 
element values are the same quantities obtained by using f? — 7. 5 
degrees. A comparison of the results entered in Tables ХІІ, 

XIII and XIV shows essentially zero error. This is well within 
normal engineering tolerances for structural analysis and dice 
that under proper conditions either of the two methods is capable 
of satisfactory accuracy. 

Since accuracy alone does not provide a sufficient basis for 
choosing one method ower thevother, the comvembeñce or eğse of 
usage of each method must be examined. Ease of usage must be 
considered within the framework of the load applied, the accuracy 
desired, the equipment available, and the assumptions made. A 
comparison of the ease of usage of the two methods within the 
framework of the load applied reveals that the Fourier series 
method has distinct advantages over the finite element method when 
solving the ring problem for a load that is applied over the surface 
of the ring in such a manner that P p is a function that is easily 


evaluated for changing m and (. There are two reasons for this. 
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шэ wi ХИ 


RESULTS FOR CASE 1 LOAD 


0. 000 





0.318 


=0. 250 


0.068 


-0. 433 
-0. 115 


0. 0683 o- 


^^ 


-0. 500 - 0.500 


-0. 182 | збе 
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TABLE XIII 


RESULTS FOR CASE 2 LOAD 


Oe - 
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TABLE XIV 


HESULTS FOR CASE 3 LOAD 


0 Fourier Series Finite Element Analytical 
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The first is that the series 15 easily evaluated when r,, is a 
simple — The second is that for distributed loads the 
finite element method required initial fixed end conditions for 
those elements containing the loads. Unless these initial end 
cónditions are available from a handbook or other source,they 
maiie calculated using energy principles or some other method. 
As can be seen in Appendix B, calculation of these initial end 
conditions may be laborious and time consuming and may re- 
quire a high degree of initial accuracy for the results to be 
useful. Even if the equations for calculating the initial end 
гован are available in a handbook, they may be їен 0 
and may also require a high A тес of initial Usos. 

A situation in which the finite element method is superior, 
in the framework of the load applied, is one where the load con- 
ditions lead to either t = 0 initial end conditions or fixed end 
conditions that are available and easily calculated. This is 
especially significant if several such loading conditions are 
to be considered. All that is required før a solution in such 
a situation, once a computer program has been written, is to 
change the program inputs to reflect the changes in{x] : 

(х2| and Га) caused by the load change and to run the program 
for each loading condition. 1 the Fomrier series method is 
used in this situation, each load change requires an integration 
to obtain E and a different series must be calculated for each 
loading condition. 


AS pointed our earlier in this chapter, both merbo 6 


capable of giving excellent accuracy under certain conditions. 


ЭР 


Mos ts Fourier series method, accuracy !5-a:iewwed only b, the 
aaber of terms included im the summation prec ess. For dis- 
placements, Tables П, П1 апа IV show that very few terms are 
me@uczred for a high degree of accuwacy especially when a distri- 
buted load is applied. Using only 4 terms in any of the summations 
ew, ihclüdinge that for the point load, gives a maximum error 
m O OE per cent. For No and Мф the s?tuation is not quite as 
faverable, but at the worst, 4 terms give a maximum error of 
Euer cent. if this degree of accuracy is satisfactory, the 
Fourier series method must be considered easier to use than the 
finite element method. 

For the finite element method, all element sizes considered 
gave essentialiy the same results for al! three load conditions. 
However, for load conditions other than the ones cunsidered here, | 
the analyst may be required to use elements with@ less than 7.5 
degrees. There appear to be two limiting factors on the minimum 
Eg ot Clements that can Ьё used. The frrst lumitigg factor is 
computer capacity. As the elements are made.smaller, more 
Wi are Created, aşü Pence larger matrioms are required. When 
dealime with full Scale @tructures, the limit of computer capacity 
fan be re#chdd prior to reachmygy the optimum elemeft size for 
accuracy. The cther limiting factor is related to computer 
capacity büt ariass from a characteristic ofthe curved element 
stiffness matrix rather than from the number of elements used. 

Toe characteristic referred to is the tendeacy for some of the 
numbers in the stıffness matrıces, particularly on the main diago- 
nal, to have a greater disparity in size as ( is decreased. Tables 


V, VIand VII show this to a marked degree. For instance, the 


m) 


ratio of the stiffness coefficient (6,6) to coefficient (1,1) in 

| Kn | fOr (3 = 30 degrees is 17/17,243, whereas the Bio of 

the same two elements is 69/18, 665, 685 for (3 = 7.5 degrees. 
eee the computer retains only a certain number of significant 
figures while manipulating these matrices, a point will be 
reached where round-off @rror in. compUter operations will 
negate any gain in accuracy obtained by using smaller elements. 
This problem is alleviated to a degree by using double-precision 
techniques, as used in this thesis; but again since double- 
precision operations require more computer Storage Space, the 
computer capacity itself places a limit on this approach. 

The equipment available to the analyst has a significant 
effect on the ease of usage of the two methods. Indeed, in the 
case of the finite element method, lack of a digital computer 
almost prohibits use of the method on the basis of effort required. 
Even if a digital computer is available, a desk calculator and a 
set of mathematical or trigonometric tables with up to ten signifi- 
cant figures is desirable if the nature of the load is such that 
fixed end conditions must be determined. On the other hand, 
while a desk calculator, and even a digital computer, may be used 
to advantage in the Fourier series method, neither are required.“ 
A slide-rule, a table of integrals, and normal trigonometric 
tables are sufficient equipment for application of the Fourier 
series method. 

Finally, the ease of usage must be considered within the 
framework of the theory employed. The assumptions made when 


applying the two methods of solution to the ring can have an effect 


a 


= 
Ja 
e 


on both She-accuracy obtsWedmDdlMhe-e Y т 


gc 


governing differential equations are derived, and the Fourier 
series solution developed through Equation (4-6a), (4-6b) and 
(4-6c) without making any assumptions regarding the extension 
of the middle surface of the ring. If the assumption is made 
during the derivation of these differential equations that the 
middle surface is inextensional, that is € $7 O when zz O, then 
E -W and the m= 0, or extensional, terms do not appear in 
the expressions for w and ün ™ ee seen in УГ I, to large 
error is introduced in w and M, by this assumption. However, 
if this inextensional form of w is used in Equation (2-6a) to 


dbeain-N., the resulting expression is missing the leading, or 


p’ 
m = 0, term, and the values for Ny calculated from this shortened 
expression are in unacceptable error. This problem may be 
circumvented by using the equilibrium equation, Equation (2-2a), 
to calculate h ч “x The procedure used is to substitute the inex- 
tensional form of My, obtained from Equation (2-6b), into 
Equation (2-2a)and integrate. The constant of integration is 
obtained from equilibrium considerations at one of the planes of 
symmetry. The resulting expression for No is identical to 


Equation (4-4b) and thus gives accurate results. On the other 


hand, in the finite element method, an error will occur in the 


* The governing equations for the modes m > 2 identically 
satisfy the inextensional assumption. Thus, these nodes are 
not affected by the assumption. 


** This is analgous to determining OQ, from the equilibrium 
equation, Equation (1-1a). 
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solution for w if the inextensional assumption is made and 
distributed loads are applied over almost all the ring surface. 
The reason for this is that a distributed load acting on the surface 
of a curved element has a greater tendency to compress or ex- 
tend the middle surface than does a concentrated load. Errors 
are caused if the energy of this middle surface compression or 
extension is neglected, asit 1syunsthevdemgveationwof thee an pue 
element stiffness matrix. Unlike the Fourier series method, 
there is no means of introducing this extensional behavior into 
the finite element method once the stiffness matrix is derived. 
The only recourse is to selecta t/a ratio and re-derive the 
curved element stiffness matrix including the axial strain energy 


term in the calculations. 
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CHAPTER VIII 


CONCLUSIONS AND RECOMMENDATIONS 


The primary advantage of the Fourier series method for 
analyzing ring structures is that it gives reasonably accurate 
answers with a minimum amount of effort using a slide rule 
and normally available integral and trigonometric tables. 

The primary advantage of the finite element method lies 
in its ability to accurately solve several different loading con- 
ditions with a minimum amount of effort. The finite element 
method has the disadvantage of requiring the calculation of 
fixed end conditions for certain types of loads. 

A recommendation for further work is the derivation of a 
stiffness matrix for a curved element that includes the effects 
of axial strain energy. Also recommended is a matrix error 
analysis on both the extensional and inextensional forms of the 
curved element stiffness matrix for various values of the 


central angle. 
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APPENDIX A 


ANALYTICAL SOLUTION 


The ring problem is solved analytically for No and м, for 
all three load conditions to obtain a standard to which the Fourier 
series and finite element results can be compared. The ana- 
lytical solution is based on the symmetric character of the loads 
applied, the principle of superposition and Castigliano's Theorem. 
The use of superposition in the solution of thin shell problems is 
discussed in detail in Reference 5. Figure A 1 (I) shows the 
entire ring with an applied load of P pounds uniformly distri- 
buted over each of the two surface areas subtended by the angle 
2ο. Figures А 1 (П), (III) and (IV) show the portions of the 
ring that are superimposed to obtain a solution. 

Considering the quadrant of the ring for 0= ў = 1/2, the 


boundary conditions plus equilibrium for Figure A 1 (II) require 


that 

(Np), = = B/2« (A -1а) 
and 

(М) = - 0 (A -1b) 
for 0< 0Ü < . Similarly 
q (Маш = Н,соѕ ф (А -2а) 
and 

(Мун = M, + aH, (cos Ò - cos) (A -2b) 


for Figure A 1 (111). 
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RING SEGMENTS USED IN ANALYTICA L SOLUTION 
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In Figure A 1 (IV) 


_ A Psin e sin ў 


(N hv 7 2 ёс (A -3a) 
and | 
s a Psine (l-sin ў) 
Сбу ΝΡ 2 ос 


are written fore = ф = 1T/2. Since equilibrium gives 


аР (1-ѕіпе‹ ) 


Μι E M, - all, cos © - x 


M, can be substituted in the expression for (Магу to obtain 


E d P (l-sin« si ) 
(Mo ту = M, a H, cos ex + aP ( 5-2 sin 0 (A -3b) 
When ф = =< 
so that Equation (A -3b) yields 
| Ἢ Pcose - 
Н, E a 


This expression is substituted into Equations (A -2a) and (A -2b), 
and the resulting equations are superimposed on Equations (A -la) 


and (A -1b) to obtain 


m P 
No = F (cose cos 0-1) - | (A -4a) 
and 
a P coso 
м, = M, t — (E% ў - cose) (A -4b) 


fer τος, 
Substituting the same expression into Equation (A -3b) yields 

м = м,+ 2 SİN (sino 2 sin $f (A -5) 
where x= 0 = TF /2. 


Castigliano's Theorem is applied in the form of 


Ola = | | 
м us О to obtain the internal цэвдэг M,. The total 
strain energy, Ub of the three segments is written as 
U, | 7  Ug,t Ugp! Ujrra ! Urrre!UrVa tUIVb 
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where subscripts a and b stand for axial and bending strain 

energy respectively. Theterms Um“ Зць” Ujija: and Ury} need 
not be considered when applying Castigliano's Theorem to obtain 
M, since these terms are either zero or are not functions of M3. 


тере оге Эл кы =) becomes 
a 


У << м: 2772 μι; 
© 
= — ado + a ad = 0 (A -6) 

Ум, pr pss 

О ος 

where the expressions for ma are obtained from Equations (A -4b) 
for the first integral and (A -5) for the second integral. Performing 
the partial differentiation and integration with the appropriate 


limits as indicated in Equation (A -6) gives 





а Р ff sin*ex 
М,-- π | 2 < = J (a 


Using this expression for M, in Equations (A - 4) апа (А -5), the 


2 


complete set of equations for № апа Mg is written. 


P 
Хуй гс (όσο (05 0 E 
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1 4 Cosa cos | 
= aP| — - ME — 
Ma | T 2x 2 x | 





его = < z and 


P sina sin $ 


ki 
Nb ax 


- An Sie Sin Ф 
Mas aP| й — 4 


2 = 


ба 2 ф AN AS. These equations are non-dimension- 
alized and evaluated for 0 2 0, 30, 60 and 90 degrees for all 


three case loads. The results are entered in Table AI. 
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TABLE A І 


ANALYTICA L SOLUTION RESULTS 


ΠΝ Case 1 Load Case 2 Load Case 3 Load 


0.0 “02287 20:238 






тэо 0. 080 





-0. 270 


0. 048 


-0. 358 


-0. 040 


-0.413 


-02095 
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APPENDIX B 


SOLUTION FOR INITIAL FORCES 


The techniques used to obtain the initial fixed end moments 
and forces required for the finite element solutions to Case 2 and 
Case 3 loads involve superposition and Castigliano's Theorem. 
Figure B 1 shows the finite element, the three loading conditions, 
the boundary forces and moments, and the notation used. The 
superposition of the forces and moments in Figure B 1 (II) and 
B 1 (111) is equivalent to the forces and moments in Figure B 1 (1). 
The sum of the rotations of the ends of the element in Figure 
B 1 (II) and (III) must be zero. Thus, 63; + Өтүү = © 
However, since a is zero under the end conditions and loading 
of Figure B 1 (II), 9j = 0. The analogous requirement for the 
horizontal displacements of the ends of the element of Figure 


B ! (II) and (III) leads to 


54 - r (B-1) 
where Ya is the horizontal component of the known radial dis- 
placement of (11), and Зил is the horizontal displacement obtained 
from an expression for the strain energy of the element in Figure 
B 1 (II) using Castigliano's Theorem. эн is positive in the 


direction towards the centerline, and is positive in the 


Sm 


direction away from the centerline. 
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Letting x be the half angle of the element and ( be the 


angle at any point along the element, the moment at 0 15 


Me = M?-0Hg cos O - Coso: ] 
And the axial force at Q is 


T = Wik. ces ф 


The strain energy, U, from both axial extension and bending 


| of the total. element 15 


ΜΓ Έα. 
υκΞ ΕΙ — Jb | : 





ade (B-2) 
? FS | 2АЕ 
“о 
Using Equation (B-2) in Castigliano's Theorem gives 
: I \) 
Uu rr > тө 
Ш ΝΗ ай 

Substituting Equation (B-2) into Equation (B-3), differentiating, 


integrating and applying limits yield 


^ 





A 15.1 E = = ose} re 
Зай эс! цн һе — же ση 
| A an 
py, An E e соёж| a UI (B-4) 
A RE. TE 
Sin ἐς 
+ | ) 
H j 
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The horizontal displacement S II is 


5 ? a Psine > 
З Zeviítt E 


where f is given by 


CASE LOAD 


-ar 





The factor f is required in order to give the proper portion of the 


total load P over each of the different segment sizes. Substi- 


tuting Equations (B-4) and (B-5) into (B-1) leads to 
i 


Б ч" 3 SY 2« 2 
S - ={sinx- cose} $e HK COS οί 
тї ШШ) х 2 А 


О, ex Sim ax a P sina 
H — + ——— |= —s— 
+ ae tE 


+ —— 
AE M| 2 
The above equation can be simplified to. give 


2 Ὅτ т» 
12.5 На = {s\nx-acose}+ Έα 


2 4 


2 + | 25$ 





2 ` Эс P sin < 
+x COS «x an7 ыг Ээ 2 — (B-6) 


Thus, Hy¡ycan be determined from Equation (B-6) for a given case 
load and a specified central angle for the element, B , Once t/a is 


chosen. The requirement to select a specific t/a ratio in order to 
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solve Equation (B-6) means that the axial strain energy is in- 
cluded in the expression. Failure to include the axial strain 
energy by making the assumption made in Reference Z2 that it can 
be neglected for curved elements with small t/a ratios leads to 
gross errors. This is because the assumption made in Refer- 
ence 2 is not valid for curved elements with distributed loads 
over large portions of the element. 

The equation for the total normal forces М? acting on the 


ends of the element shown in Figure B 1 (I) is 


N° = Nac Hu 695 = 


or 
o р 
NM 2 —- Ha cose (8-7) 
2e 
since p 
№ = 2 — £ 
The component of H,., in the radial direction is the shear 


IH 


force, 09, at the ends of the ring element. Thus 


GS Sec (B-8) 


In order to obtain M?, Castigliano's Theorem, in con- 


junction with the requirement 6,,,; = 0, is used to write 


ЭМ o 
3M? 


where U, is given by Equation (B-2). 
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Differentiation and integration of this equation give 


SMS. -- E 
О 


Experience with ps (B-6) through (B-9) shows that the 
axial strain energy must be included if accurate answers are de- 
sired [0 №, Q% and MP. Thus, t/a 2 1/10 is chosen, and 
Equations (B-6) through (B-9) are used to prepare the entires 
in Table B 1 for NF i and M? for Case 2 and 3 loads and f? = 30, 
15 and 7. 5 RN In evaluating these equations, trigonometric 
functions accurate to 5 decimal places do not give accurate 
answers for N?,. Q? and M? for the values of (3 used in this 
thesis. Consequently the entries in Table B Iare obtained by 
using trigonometric functions of 10 decimal place accuracy. 
Solving equations such as (B-6) --- decimal places is ob- 
viously beyond the capacity of either slide rule or practical hand 


calculations. Even with an electronic desk calculator, solving 


such equations is laborious and time consuming. 
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TABLE B I 
FIXED END RESULTS 


CASE 2 LOAD 


NS --——— - 
| тээ 15 


nd 


0. 00181929P 0. 01286764P 








| 0. 12707678Р 
9” 0. 06247015P 0. 12402482P | 0. 22182251Р 
| 


0. 00136328Ра 0. 0054178 0Ра 0. 01944670Ра 


CASE 3 LOAD 


132 44 | 30° 


——— ---- 








0. 00090964P 0. 00643382P 0. 06353839P 






"0. 03123508Р 0. 06201237P 0. 11091125P 






0. 00068 164Pa 0. 0027089 0Pa 





0. 00972335 Pa- 
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А ЕРЕКЕ 
FORTRAN IV PROGRAMS 
TABLE C I 


DESCRIPTION OF SUBROUTINE "DSIMQ" 


PURPOSE 
OBTAIN SOLUTION OF A SET OF SIMULTANEOUS LINEAR 
EQUATIONS, AX=B 


USAGE 
CALL "DSIMQ" (A,B,N, KS) 


DESCRIPTION OF PARAMETERS 

A AND B MUST BE REAL*8 

A - MATRIX OF COEFFICIENTS STORED COLUMNWISE. 
THESE ARE DESTROYED IN THE COMPUTATION. 
THE SIZE OF MATRIX A IS N BY N. 

B.- VECTORS OF ORIGINA L CONSTANTS (LENGTH N). 
THESE AREREPLACED BY FINAL SOLUTION VALUES, 
VECTOR X. 

N - NUMBER OF EQUATIONS AND VARIABLES 

KS -OUTPUT DIGIT 

| 0 FOR A NORMAL SOLUTION 
] FOR A SINGULAR SET OF EQUA TIONS 


REMARKS 
MATRIX A MUST BE GENERAL, 
IF MATRIX IS SINGULAR, SOLUTION VALUES ARE MEAN- 
INGLESS. AN ALTERNATIVE SOLUTION MA Y BE OB- 
TAINED BY USING MA TRIX INV ERSION (MINV ) AND 
MATRIX PRODUCT (GMPRD). 


SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED 
NONE 


METHOD 
METHOD OF SOLUTION IS BY ELIMINATION USING 
LARGEST PIVOTAL DIVISOR. EACH STAGE OF ELIMINA- 
TION CONSISTS OF INTERCHANGING ROWS WHEN NECES- 
SARY TO AVOID DIVISION BY ZERO OR SMALL ELEMENTS. 
THE FORWARD SOLUTION TO OBTAIN VARIABLE N IS 
DONE IN N STAGES. THE BACK SOLUTION FOR THE 
OTHER VARIABLES IS CALCULATED BY SUCCESSIVE 
SUBSTITUTIONS. FINAL SOLUTION VALUES ARE 
DEVELOPED IN VECTOR B, WITH VARIABLE 1 IN B(1), 
VARIABLE 2.IN B(2) . . . . . . , VARIABLE N IN B(N). 
IF NO PIVOT CAN BE FOUND EXCEEDING A TOLERANCE 
OF 0.0, THE MATRIX IS CONSIDER ED SINGULAR AND 
KS IS SET TO 1. THIS TOLERANCE CAN BE MODIFIED 
BY REPLACING THE FIRST STATEMENT. 


ME 


TABLE C II 


LIST OF SYMBOLS USED IN SUBROUTINE "RELM" 


Computer 
Coded Name 


DBETA 
DESM 


DESM (I, J) 


DLA 
DLB 
DLC 
DLD 
DLE 
DBA 
DBB 
DBC 
DBD 
DBE 
DBF 


DBG 


Definition 
/5 , element central angle. 
Element stiffness matrix [К η! : 


Stiffness influence coefficient in 
K 
[Ка 


а (corresponding notation in 
b Reference 3). 
с 

d 

e 

A 

B 

C 

D 

E 

P 

G 


73 


5 


TABLE CIII 
SUBROUTINE "RELM" 
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TABLE C IV 
LIST OF SYMBOLS USED IN 


PROGRAM "BETA 30" 


Computer 
Coded Name Definition 
ы Ὁ 
А Reduced stiffness matrix Kp. 
A (I, J) Stiffness influence coefficient 
in | Koj. 
B (I) Prior to calling SDSIMQ'' B(I) = 
v oa ον) 
B (1) After calling "DSIMQ" B (I)= {up}. 
DBETA (4 , element central angle. 
DESM Element stiffness matrix b.) 
DESM (I, J) Stiffness influence coefficient 
in ES . 
L 
PI fT 
DM (I) Moment at node (I). 
DN (I) Normal force at node (I). 
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required. The results indicate a satisfactory accuracy for both 
methods under most conditions. The Fourier series method is 

! superior for solving problems with a distributed load condition. The 
finite element method is superior for solving multiple problems with 
concentrated loads. 
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